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NONSYMMETRIC  BALLASTIC  RANGE,  HEIGHT,  TIME-OF-FLIGHT 
AND  OPTIMAL  FLIGHT  PATH  ANGLE  COMPUTATIONS 
WITH  PROGRAMS  FOR  A  HEWLETT-PACKARD  65  CALCULATOR 


I .  The  Introduction  and  Apology. 

The  range  of  a  ballistic  missile  over  the  surface  of  the 
Earth  is  a  topic  which  is  discussed  in  almost  every  text  on 
astronautics  (Reference  1,  for  example).  Unfortunately,  none 
of  the  available  texts  contained  a  needed  procedure  for  the 
exact  solution  of  the  unsymmetric  ballistic  problem  (defined 
in  Section  II) ;  all  procedures  found  were  approximate  and 
were  based  on  a  symmetric  approximation  [See  References  2  and 
3/  for  example] .  The  author  of  this  report  faced  the  dilemma 
of. either  (1)  taking  the  time  to  solve  the  problem  afresh,  or 
(2)  taking  an  inordinate  longer  time  to  do  a  thorough  search 
to  find  a  solution  which  might  not  be  the  one  desired.  The 
iormer  course  of  action  was  taken  and  all  due  apologies  are 
hereby  extended  to  those  that  have  published  more  elegant 
solutions . 

All  of  the  computational  methods  are  outlined  in  Section 
III.  Instructions  and  a  sample  problem  for  the  Hewlett- 
Packard  65  Calculator  (HP-65)  are  given  in  Section  IV.  This 
HP-65  programs  are  in  Section  V,  and  the  development  of  the 
3opt  procedure  is  in  the  Appendix.  Familiarity  with 
Newtonian  two-body  theory  is  assumed. 


II.  The  Problem 


Basically,  the  problem  is  to  compute  the  range,  S,  and 
time-of-f light ,  t^,  of  a  ballistic  missile  given  the  booster 
cutoff  at  height  h„  after  a  vertical  ascent,  a  cutoff 
velocity  v^^,  and  a  subsequent  flight  path  angle  6  (Fig.  1)  . 
Given  that  a  solution  to  the  stated  problem  is  available,  the 
next  natural  question  to  ask  is,  "What  flight  path  angle 
will  obtain  the  maximum  range  S  ?"  The  solution  is  fairly 
simple  in  the  symmetric  ballistic  problem,  i.e,  when  h^  =  0, 
and  called  symmetric  because  the  launch  point  L  and  target 
point  T  are  symmetric  with  the  line  of  apses  of  the 
Keplerian  orbit  (i.e.  (|)j^  =  (f)2  in  Fig.  1).  The  unsymmetric 

ballistic  problem  is  that  in  which  h_  ^  0.  Here  we  assume 

£j 

that  hg  is  given,  the  Earth  with  radius  r^  is  spherical, 
non-rotating,  and  has  no  atmosphere.  The  effect  of  the  Earth's 
rotation  may  be  calculated  figuring  how  far  the  target  has 
moved  during  the  time  of  flight  of  the  missile.  Explicitly, 
we  assiame  that  Keplerian  orbit  is  obtained  at  height  h_  with 

rj 

terminal  booster  cutoff  velocity  v_  and  flight  path  angle 
6  (or  ,  and  the  target  T  is  at  the  Earth's  surface. 

The  computation  of  range  and  time-of-f light  for  given  h^, 
^E'  ^  easily  be  obtained  from  the  equations  of 

elliptical  orbits.  The  determination  of  is  somewhat 

complicated,  however,  and  requires  an  iterative  procedure.  In 
all  of  the  developments  we  assume  that  L  and  T  are  sepa¬ 
rated  by  the  line  of  apses  (which  implies  that  0  <  &  <  ^  ) . 
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Line  of  apses 


Unpowered  flight 


FIGURE  1 


III.  The  Computational  Procedures. 


A.  Range 


1.  Let  and  denote  the  distance 

of  the  launch  point  L  and  target  point 
T  from  the  center  of  the  Earth.  Usually 
we  let  rj^  =  r^+h^  and  r^  =  r^ 
although  and  r^  need  not  be  so 

defined  in  general. 

2 .  From  r^  and  v,.  ,  compute  the  semimajor 

Li  Li 

axis  of  the  ellipse. 


a 


2 

where  y  =  and  is  the  (Gaussian) 


gravitational  constant  for  the  Earth.  We 
have  assumed  that  the  mass  of  the  missile 
is  negligible  with  respect  to  the  mass  of 
•  the  Earth. 

3.  Using  the  flight  path  angle  6,  compute 
the  elliptical  semiparameter 


4.  Compute  the  eccentricity 
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5.  Compute  the  true  anomaly  of  the  launch 


point 


where  0  <  f_  <  ir  .  Note  that  the  line 

Li 

of  apses  can  be  said  to  have  true  anomaly 

f  =  IT  . 

6 .  Compute  the  true  anomaly  of  the  target  point 


7.  The  surface  range  of  the  missile  is  then 

S  =  r^(f^  -  fj.)  . 


B.  Time-of -Flight 

1.  Compute  the  eccentric  anomaly  of  the 
launch  point. 


where  0  <  E  <  ir  . 

Li 

2 .  Compute  the  eccentric  anomaly  of  the 
target  point. 


where  it  <  <  27r  . 
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3.  Using  Keplers'  equation,  compute  the 


time-of -flight 


tf  = 


/vT 


(E^  -  Ej^)-e{sinE^  -  sinEj^) 


C. 


Height  of  Trajectory  &  Circular  Velocity 

1.  The  maximum  height  of  the  trajectory  above 
the  Earth's  surface  for  any  h^,  and 

3  is  given  by 


Nax  =  h  =  a(l  +  e)  - 


2.  Circular  Velocity 

The  velocity  required  for  a  circular 

orbit  at  r,.  is 
Jb 


’'c  = 


V  will  be  used  in  Section  III  D, 
c 


D.  Optimum  Flight  Path  Angle  &  Maximum  Range 

For  a  specified  h„  and  v,.  ,  a  value  of 

£i  Li 

3  can  be  found  such  that 


Smax  =  S'Vt>  '  • 


In  this  unconstrained  optimization  problem 


we  set 


93 


3  =  3 


=  0 


and  solve 


opt 
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for  relationship  between  S  and 

B  is  given  by  the  computational  formulas  for 
the  range  in  part  A.  Although  an  explicit 
formula  for  $  is  unobtainable,  the  following 
exact  relationship  is  useful: 


where  a  =  —  . 


In  the  symmetric  ballistic  problem  r^^  =  r^  ; 
the  equation  for  ^Qp^.  then  reduces  to 


(2) 


cos^B 


opt 


2-a 


2  ' 


which  is  the  result  given  in  References  1,  2,  and  3. 


find 


In  the  unsymmetric  case  Equation  (1)  is  used  to 

B^  .  as  follows: 
opt 


1.  Solve  (2)  for  B  • 

opt 

2.  Substitute  B  into  the  right  hand  side 

opt  ^ 

of  (1)  to  obtain  a  new  value  for  B  a.  • 

opt 
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3.  Average  the  last  two  values  of  to 

obtain  a  new  value  of  and  return 

opt 

this  value  of  to  step  2. 

4.  It  has  been  found  that  five  iterations  are 
sufficient  to  obtain  a  reliable  value  for 

^opt 

5.  Using  conipute  the  maximum  range 

^max  “  ^^^opt^  using  procedure  A. 
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IV.  The  Illustrative  Example. 


The  HP-65  program  consists  of  four  magnetic  cards  which 
are  used  as  follows: 

Card  1:  This  card  must  be  entered  first.  It  is  used 
to  set  the  two  constants  r„  and  y  for  any 

III 

choice  of  four  distance  units  and  any  choice  of 
four  time  units.  The  distance  units  are 
kilometers,  meters,  statute  miles,  and  nautical 
miles;  the  distance  unit  is  entered  by  pressing 
key  A,  B,  C,  or  D,  respectively.  The  display 
shows  the  equatorial  Earth  radius  r^  in  the 
respective  unit.  The  time  units  are  seconds, 
minutes,  hours,  and  days;  the  time  unit  is 
entered  by  pressing  key  A,  B,  C,  or  D,  respec¬ 
tively.  The  display  shows  the  constant  y  in 

3.2. 

the  appropriate  (distance)  /(time)  unit.  Note  - 
after  entering  program  Card  1,  the  distance  unit 
must  be  entered  before  the  time  unit  is  entered 
since  the  keys  A  through  D  perform  a  dual  func¬ 
tion.  The  basic  definitions  used  are: 

1  statute  mile  =  1.609344  kilometers (km) 

1  nautical  mile  =  1.852  km 
r^  =  6378.160  km 

Ili 

y  =  398603  km^/sec^ 

Card  2:  This  card  computes  S,  h  ,  and  v  given 

lUaX  O 

h„,  and  3.  It  must  be  used  after  the 
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execution  of  Card  1  and  it  may  be  used  before 
and/or  after  the  execution  of  Cards  3  and  4.  If 
it  is  used  after  Card  4  to  compute  then 

as  computed  by  Card  4  need  not  be  re-entered. 
Otherwise,  in  any  order:  enter  h„  and  press  A, 

Ij 

enter  and  press  B,  enter  6  and  press  C.  To 
find  the  range  S,  press  E.  Then  to  find  ^max' 
press  R/S;  and  then  to  find  v^,  press  R/S  again. 
Note;  If  h_,  Vt ,  and  3  have  been  entered 
using  Card  3 ,  then  they  do  not  have  to  be 
reentered  when  using  Card  2.  Also,  any  new  value 
of  hg,  Vj^,  or  3  can  be  entered  and  the  range 
computed  without  having  to  re-enter  unchanged 
values. 

Card  3:  This  card  computes  the  time  of  flight,  t^  . 

It  must  be  used  after  the  execution  of  Card  1 
and  it  may  be  used  before  and/or  after  the 
execution  of  Cards  2  and  4.  If  it  is  used 
after  Card  4 ,  then  computed  by 

Card  4  need  not  be  re-entered;  if  it  is  used 
after  Card  2,  then  only  new  quantities  need  be 
entered.  Otherwise,  in  any  order:  enter  hg 
and  press  A,  enter  Vg  and  press  B,  enter  3 
and  press  C.  To  find  the  time  of  flight,  press 
E. 
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Card  4;  This  card  computes  ^Qp^*  must  be  used 

after  and  have  been  entered  from 

either  Card  2  or  Card  3.  No  entries  are 
required.  To  compute  P^^ess  any  key  A 

through  E;  five  iterations  will  be  performed 
and  the  resulting  value  of  ^Qp^  will  be  dis¬ 
played.  To  perform  five  more  iterations,  press 
R/S.  To  monitor  the  intermediate  iterations 
press  keys  £  SF2 .  Each  intermediate  result  will 
be  displayed.  To  continue,  press  R/S.  To  dis¬ 
able  the  monitoring,  press  keys  f  ^  SF2 . 

Using  the  four  program  cards  as  described  above  find 
'  the  range  and  time-of-flight  for  a  ballistic  missile  that 
has  a  velocity  of  10000  knots  at  a  cut-off  altitude  of  50 
nautical  miles  for  a  flight  path  angle  of  30*^  and  40*^. 
Also,  find  the  maximum  possible  range  and  time-of-flight 
for  the  same  missile. 


Proceed  as  follows: 


Instruction 

Data 

Key 

Display 

Enter  Card  1 

Distance  unit  is  n.mi. 

D 

rg=3443.93n.mi 

Time  unit  is  hours 

C 

y=8. 13247x10^^ 

Enter  Card  2 

Enter  cutoff  altitude 

50 

A 

r.  =3493 . 93n .mi 

Li 

Enter  cutoff  velocity 

10000 

B 

Enter  flight  path  angle 

30 

C 

Compute  range 

E 

S=1928 . 08n.mi . 
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Step 

Instruction 

Data 

Key 

Optional:  Display 

luaX 

R/S 

8b 

Optional:  Display  v 

c 

R/S 

9 

Enter  flight  path  angle 

40 

C 

10 

Compute  range 

E 

11 

Enter  Card  3 

12 

Compute  time-of-flight 

E 

Note:  hg  ~  50,  =  10000,  and 

3  =  40°  since  no  further 

entries  have  been  made. 

13 

Enter  flight  path  angle 

30 

C 

14 

Compute  time-of-flight 

E 

15 

Enter  Card  4 

16 

Compute  ^opt  iterations) 

E 

17 

Make  5  more  iterations  to  be  sure 

R/S 

18 

Enter  Card  2 

19 

Compute 

max 

E 

20 

Enter  Card  3 

21 

Compute  time-of-flight  for 

E 

maximum  range. 

Display 

=344.74n.mi 

max 

V  =15256. 47knot 
c 

S=1956.60n.roi. 
t^=0 . 3043hrs, 


0 . 2504hr3 . 


‘f<8opt>">-28« 


} 
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V.  The  Program  Listings. 
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CARD  1 


CARD  2 


Set  units  of  distance  and  time 
display  r_  and  y 


Compute  range , 


10 


20 


30 


40 


Key 

Code 

Key 

Code 

Key 

Cpde 

Key 

Code 

Entry 

Shown 

Entry 

Shown 

Entry 

Shown 

Entry 

Shown 

LBL 

23 

3 

03 

RCL  4 

3404 

X 

71 

B 

12 

7 

07 

+ 

61 

1 

01 

EEX 

43 

8 

08 

STO  3 

3303 

- 

51 

3 

03 

• 

83 

R/S 

84 

CHS 

42 

CHS 

42 

1 

01 

LBL 

23 

f 

31 

GTO 

22 

6 

06 

B 

12 

y - 

09 

0 

00 

X 

71 

STO  2 

3302 

STO  6 

3306 

LBL 

23 

STO  4 

3304 

R/S 

84 

RCL  8 

3408 

C 

13 

R/S 

84 

LBL 

23 

RCL  3 

3403 

1 

01 

60 

LBL 

23 

10 

C 

13 

60 

D 

14 

• 

83 

B 

12 

STO  7 

3307 

RCL  8 

3408 

6 

06 

6 

06 

R/S 

84 

RCL  4 

3404 

0 

00 

0 

00 

LBL 

23 

D 

14 

9 

09 

GTO 

22 

D 

14 

+ 

61 

3 

03 

1 

01 

g 

35 

g 

35 

4 

04 

LBL 

23 

RAD 

42 

IT 

02 

4 

04 

C 

13 

81 

ENTER 

41 

GTO 

22 

3 

03 

1 

01 

+ 

61 

0 

00 

6 

06 

- 

51 

g  x<— >y 

3507 

LBL 

23 

70 

0 

00 

20 

RCL  6 

3406 

70 

51 

D 

14 

0 

00 

81 

RCL  4 

3404 

1 

01 

GTO 

22 

f-1 

32 

X 

71 

• 

83 

1 

01 

COS 

05 

R/S 

84 

8 

08 

LBL 

23 

RTN 

24 

RCL  6 

3406 

5 

05 

D 

14 

LBL 

23 

1 

01 

2 

02 

•  8 

08 

E 

15 

+ 

61 

GTO 

22 

6 

06 

g 

35 

RCL  5 

3405 

0 

00 

4 

04 

DEG 

41 

81 

LBL 

23 

0 

00 

RCL  7 

3407 

RCL  4 

3404 

A 

11 

80 

0 

00 

30 

f 

31 

80 

51 

1 

01 

GTO 

22 

COS 

05 

R/S 

84 

LBL 

23 

1 

01 

RCL  2 

3402 

RCL  1 

3401 

0 

00 

LBL 

23 

X 

71 

RCL  3 

3403 

g 

35 

A 

11 

RCL  3 

3403 

i 

81 

1/x 

04 

1 

01 

X 

71 

f 

31 

ENTER 

41 

LBL 

23 

ENTER 

41 

- - 

09 

ENTER 

41 

1 

01 

ENTER 

41 

ENTER 

41 

X 

71 

R/S 

84 

X 

71 

X 

71 

RCL  1 

3401 

GTO 

22 

X 

71 

90 

RCL  1 

3401 

81 

E 

15 

3 

03 

X 

71 

40 

STO  8 

3308 

90 

g  NOP 

3501 

9 

09 

STO  1 

3301 

2 

02 

g  NOP 

3501 

8 

08 

R/S 

84 

RCL  3 

3403 

II 

3501 

6 

06 

g  NOP 

3501 

T 

81 

II 

3501 

0 

00 

g  NOP 

3501 

RCL  2 

3402 

tl 

3501 

3 

03 

II 

3501 

ENTER 

41 

It 

3501 

X 

71 

II 

3501 

X 

71 

II 

3501 

STO  1 

3301 

II 

3501 

RCL  1 

3401 

n 

3501 

CL  X 

44 

II 

3501 

T 

81 

II 

3501 

6 

06 

100 

II 

3501 

- 

51 

II 

3501 

1 

1 

50 

STO  5 

3305 

IOC 

II 

3501 

NOTE:  Before  putting  a  program  into  memory,  after  switching  to  W/PGM  mode,  press  f 
and  then  PRGM  to  clear  memory. 


CARD  3 


CARD  4 


Compute  time-of-f light 


Key 

Entry 

STO  5 
X 
1 

CHS 

f 


Key 

Code 

Entry 

Shown 

RCL  4 

3404 

+ 

61 

STO  3 

3303 

R/S 

84 

LBL 

23 

B 

12 

STO  2 

3302 

R/S 

84 

LBL 

23 

C 

13 

STO  7 

3307 

R/S 

84 

LBL 

23 

D 

14 

g 

35 

RAD 

42 

RCL  5 

3405 

X 

71 

1 

01 

- 

51 

CHS 

42 

RCL  6 

3406 

-1 

81 

f 

32 

COS 

05 

RTN 

24 

LBL 

23 

E 

15 

g 

35 

DEG 

41 

RCL  7 

3407 

f 

31 

COS 

05 

RCL  2 

3402 

X 

71 

RCL  3 

3403 

X 

71 

ENTER 

41 

X 

71 

RCL  1 

3401 

A 

81 

STO  8 

3308 

2 

02 

RCL  3 

3403 

81 

RCL  2 

3402 

ENTER 

41 

X 

71 

RCL  1 

3401 

1 

81 

60 


70 


80 


90 


/ 

STO  6 
g 

TT 

ENTER 

+ 

RCL  4 
D 

ENTER 

f 

SIN 
RCL  6 
X 


RCL  3 
D 

ENTER 

f 

SIN 
RCL  6 
X 


RCL  5 

RCL  1 
RCL  5 
X 
f 


R/S 

CTO 

E 

g  NOP 
g  NOP 
g  NOP 


100 


Code 

Shown 

51 

3305 
71 
01 
51 
42 
31 
09 

3306 
35 
02 
41 
61 

3404 
14 
51 
41 
31 
04 

3406 

71 

51 

3403 

14 
41 
31 
04 

3406 

71 

51 

51 

3405 
81 

3401 

3405 

71 

31 

09 

81 

84 

22 

15 
3501 
3501 
3501 
3501 
3501 
3501 
3501 
3501 


Compute  maximizing  flight  path  angle  8 


10 


20 


30 


40 


50 


Key 

Entry 

2 

RCL  2 
ENTER 
X 

STO  5 
RCL  3 
X 

RCL  1 


STO  6 

RCL  4 
g 

1/x 

RCL  3 
g 

1/x 

RCL  5 

RCL  1 
X 

ENTER 

+ 

1 

+ 

STO  5 
g 

DEG 

0 

STO  7 
5 

STO  8 

LBL 

0 

RCL  5 

RCL  3 

RCL  4 

RCL  6 
X 

RCL  7 
f 

COS 

X 

gLSTx 

X 

RCL  7 
f 


Code 

Shown 

02 

3402 
41 
71 

3305 

3403 
71 

3401 

81 

51 

3306 

3404 
35 
04 

3403 

35 

04 

51 

3405 
81 

3401 

71 

41 

61 

01 

61 

3305 

35 

41 

00 

3307 
05 

3308 
23 
00 

3405 

3403 

3404 
81 

3406 
71 

3407 
31 
05 
71 

3500 

71 

51 

3407 

31 


60 


70 


80 


90 


100 


Key 

Entry 

SIN 

X 

RCL  5 
RCL  3 
RCL  4 

RCL  7 
f 

COS 

X 

ENTER 

X 


1 

+ 

RCL  6 


COS 

f 

TF  2 
R/S 
g  NOP 
RCL  7 
gx<— >y 
STO  7 
0 

g  x=y 
GTO 
0 
+ 

+ 

2 

STO  7 
g 

DSZ 

GTO 

0 

R/S 

5 

STO  8 
GTO 
0 

g  NOP 


opt 

Code 

Shown 

04 

71 

3405 

3403 

3404 
81 

3407 

31 

05 

71 

41 

71 

51 

31 

09 

81 

01 

61 

3406 
81 

31 
09 

32 

05 

31 

81 

84 

3501 

3407 
3507 
3307 

00 

3523 

22 

00 

61 

61 

02 

81 

3307 
35 

83 
22 
00 

84 
05 

3308 
22 
00 

3501 


NOTE:  Before  putting  a  program  into  memory,  after  switching  to  W/PGM  mode,  press  f 

and  then  PRGM  to  clear  memory . 
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VII.  The  Appendix. 


As  indicated  in  section  HID,  the  maximum  range  for  a 

specified  cutoff  altitude  h„  and  cutoff  velocity  is 

h  L 

found  by  the  following  unconstrained  optimization: 


0<6s 


IT 


Using  the  equations  in  section  IIIA,  we  want  to  find 

3  .  such  that 

opt 


as 

aB 


aB 


B  =  B 


=  0 


opt 


B  =  B 


opt 


We  find  that 


aB 


where 


7  W- 


II 

1 

^  , 

aB 

1  +  1 

/p_  _  i\ 

er .  o  3 

1  2ae 

l^i  ) 

respect 

,  i=T,L 


vanish  at  B  =  B  . /  we  must  have  either: 

opt 


1. 


B  =  B  =0 
opt 


15 


or 


2. 


iP 

ae 


0 


6  =  3 


opt 


We  find  that 


2  2 

r  sin 

i-J  -U _ 

y 


23 


0 


TT 

for  3  =  0,  ±  /2,  ±  TT,  etc.  Since  we  require  that  the 

launch  point  and  the  target  point  be  separated  by  the  line 

TT 

of  apses,  we  require  that  e<3^-2  •  This  requirement  is 
imposed  because  we  must  properly  interface  the  disconti¬ 
nuity  caused  by  the  principal  angle  in  the  arc  cosine  used 
in  steps  5  and  6  of  section  III  A.  Thus  3  =  ^/2  is  the 
only  angle  of  interest  for  which  3p/33  =  0,  and  this 
angle  clearly  gives  rise  to  a  minimum  in  range  since  the 
trajectory  of  the  missile  is  vertical  (straight  up,  and 
then  straight  down) . 

The  value  of  3  which  maximizes  S  is  found  from  the 
requirement  that 


(Ym-YT ) a  -  a  =0 

^  ^  ^  -  ^opt 


This  equation  has  defied  attempts  to  find  an  analytic 
solution,  but  it  can  be  rearranged  and  displayed  as 
equation  1  in  section  III  D,  where  the  solution  procedure 
is  detailed. 


16 


INITIAL  DISTRIBUTION  LIST 


Dean  of  Research,  Code  023 
Naval  Postgraduate  School 
Monterey,  California  93940 


Defense  Documentation  Center 
Cameron  Station 
Alexandria,  Virginia  22314 


Library,  Code  0212 

Naval  Postgraduate  School 

Monterey,  California  93940 


Library,  Code  55 

Department  of  Operations  Research 
and  Administrative  Sciences 
Naval  Postgraduate  School 
Monterey,  California  93940 


Code  55Su 

Naval  Postgraduate  School 
Monterey,  California  93940 


